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Abstract 



We calculate the photon production cross-section in pA collisions un- 
der the assumption that the nucleus has reached the saturation regime, 
while the proton can be described by the standard parton distribution 
£L|| functions. We show that due to the strong classical held 0(l/g) of the 

nucleus, bremsstrahlung diagrams become dominant over the direct pho- 
' ton diagrams. In particular, we show that 7— jet transverse momentum 

, spectrum and correlations are very sensitive to gluon saturation effects in 

the nucleus. 
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1 Introduction 



An outstanding question in the description of hadronic interactions is the prob- 
lem of the parton distribution functions at very small values of the parton mo- 
mentum fraction x (see |l], ^ for a pedagogical introduction). Indeed, it is well 
known that the solution of the (linear) BFKL equation [|, extended to 

very small values of x, leads to cross-sections that increase asymptotically like 
a power of the center of mass energy, in contradiction with bounds coming from 
unitarity considerations. It is expected that saturation effects must limit the 
growth of parton distributions at small values of a; @, ||] ■ It has been argued 
that this description becomes inadequate since non-linear effects should become 
important at small x when the corresponding occupation numbers reach a value 
comparable to l/a s ||, [l(j, [Tl|] . 
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A theoretical description of this phenomenon is provided by the Colored 
Glass Condensate model jl2, 13, 14J. In this model, small- a; gluons are described 
by a classical color field, due to the fact that these modes have a large occupation 
number. This color field is driven by a classical Yang-Mills equation whose 
source term is provided by faster partons. As one varies the separation scale 
between the soft and hard modes, one obtains a renormalization group equation 
that provides a non-linear generalization of the BFKL evolution equation p5| , 
[jj] |J, H M 111 HI' if H H HI IZl> !§■ An important parameter of the 
model is the saturation momentum scale, usually denoted Q s and defined as the 
transverse momentum scale below which saturation effects become important. 
This scale increases with energy and with the size of the nucleus @, |, 0. 

More recently, some tests have been proposed in order to observe saturation 
effects in high energy hadronic or nuclear collisions. Among them are the re- 
analysis of HERA deep inelastic results under the hypothesis that the 7*-p 
cross-section is affected by saturation effects [$(], [31 . In the field of heavy 
ion collisions, it has been suggested that one could predict the multiplicity by 
describing the two nuclei with classical color fields |32|, |33|, |34|, |35|, This 
idea lead to a successful description of the observed multiplicities at RHIC 
[ p?| [j8| , |39| which is however still under debate Q. It has also been proposed 
recently that saturation effects could be observed in the cleaner environment 
provided by ultra-peripheral heavy ion collisions El], [l2| since the production 
of heavy quark pairs is also very sensitive to the saturation scale. 

Various studies also proposed to study saturation effects in the context of 
pA collisions |43, 44 45|. In this paper, we calculate the cross-section for photon 
production in pA collisions and show that the result is also sensitive to saturation 
effects. We assume that the nucleus is large enough so that it can be described 
in terms of classical color fields. Concerning the proton, we do not make such 
an assumption and describe it as in ordinary perturbative calculations, in terms 
of parton distribution functions. In this framework, the leading order process 
for photon production involves a quark from the proton that scatters off the 
color field of the nucleus and emits a photon by bremsstrahlung. Therefore, the 
elementary process we consider in this paper is q — > 57 in the background of a 
classical color field. 

The paper is organized as follows. In section ^, we formulate the amplitude 
for the bremsstrahlung of a photon by a quark (or antiquark) in the presence of 
a classical color field. This amplitude is then evaluated. The average over the 
configurations of the hard sources is performed in section ||. In section f|, we 
explain how to obtain the differential cross-section from the semi-classical am- 
plitude derived in section ||. In section ||, we show that the diffractive (i.e. with 
the additional constraint that the exchange between the quark and the nucleus 
is color singlet and carries a zero total transverse momentum) cross-section for 
this process is vanishing. Finally, in section Q we derive the inclusive photon 
radiation differential cross-section and show that the correlations between the 
outgoing quark and the photon are affected by saturation effects. 
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2 Bremsstrahlung amplitude 

We want to calculate the following processes in pA collisions 

q(p) + A^q(q) 1 (k)X (1) 
q(p) + A^q(q) 1 (k)A (2) 

where the quark can become either a jet or a hadron. 
The starting point is the amplitude: 

(q(q)j(k) out \q(p)in) = (0 O ut|aout(fc)&out(g)&L(P)|°i») ( 3 ) 
which, using the LSZ formalism Jl6| can be written as 

<0out|aout(fc)6 O ut(q)6L(p)|0in) - / ' d*xd*yd*ze i V e - x+ ''-*-™) 

^2V ^3 J 

xu(q)(i Z -m)(0 ou t\Tip(z)s- J(x)ip(y)\O in )(i y +m)u(p) 

(4) 

where J^(x) = ^}(x)^^{x) is the quark component of the electromagnetic cur- 
rent, Z2 and Z-j, are the fermion and photon wave function renormalization 
factors^], and e is the fractional electric charge of the quark under consideration. 
So what we basically need is 

(0 out | T^(2#(a;)^(a#(y) 1 in ) (5) 

which, using Wick's theorem, can be written as 

(0 OVLt \T^(z)^(x)^(x)^{y)\O in ) = ~G F {z,x)^G F {x,y) (6) 

where the fermion Feynman propagator G F is defined as[] 

G P (x,y) = <0 out |T^(y)^(x)|O in ) . (7) 

It should be noted that G F is the fermion propagator in the background of the 
strong classical color field of the nucleus. Our amplitude then becomes 



<<z(g)7(*0out|<z(p)i„) =e / d 4 xd i yd 4 ze i(k - x+ «- z -r-y'> 

Xu(q)(i z -m)G F (z,x)^G F (x,y)(i V -m)u(p) . (8) 



1 Since we are not including loop corrections in this calculation, we will simply set these 
factors to 1 in the following. 

2 Compared to the definition of the fermion propagator in ^7j, we are omitting the 
denominator (0 O ut|0j n ) due to the fact that this factor is a pure phase in the case where the 
background contains the classical field of a single nucleus. This phase indeed drops out of 
physical quantities obtained after squaring amplitudes. 
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The fermion propagator G F is already known from ^8], |49). In co- 

ordinate space and in the "singular" gauge, it is given by the following formula 



G F (x,y) = G° F (x-y) + d 4 z6(z-) 



-e(-x-)6(y-)(U(z ± )-l) 



G (x-zh-G°(z~y) , (9) 



where G° F is the free Feynman propagator of a quark, and where U(z±) is a 
unitary matrix containing the interactions between the quark and the colored 
glass condensate. It is convenient to separate out the interaction part of the 
propagator and write it in momentum space as 



G F {q,p) = (27rrsHq-p)G F (p)+Gl( q )T F ( q , P )G° F (p) 



(10) 



where 



T F (q,P) = 27r<% -p h sign(p )Jd 2 z ± 



jjsign(p )/. 



e i(q±-p±)-z ± ^-q) 



while for a nucleus moving in the positive z direction we have 



U(x±) = Texp 



dz 



1 

vi 



p a (z ,z±)t a 



(12) 



with t a in the fundamental representation, and where p a (z~ , z±) is the density 
of color sources in the nucleus. 

Inserting this propagator in the expression for the photon production am- 
plitude (H) , we find the following expression in momentum space 



(g(<zh(fc)out|g(p)in) = -eu(q) 



(27r) 4 <5 4 (fc + g- 



+T F (q,p-k) (%(p-k)j + 4Gl(q + k)T P (q + k,p) 
d 4 l 



(27T)- 



-T F (q, I) G° F (0 iG% (k + I) T F (k + l,p) 



u(p) 



(13) 



The first term in ( [13] ) corresponds to emission of a photon by the quark with 
no scattering from the Color Glass Condensate. This term cannot contribute 
to the cross section for real photons and indeed vanishes (the delta function 
8{p — q — k) has no support for on-shell particles). 

The second and third term, respectively, correspond to the case when the 
quark multiply scatters from the Color Glass Condensate before and after emis- 
sion of the photon. The last term describes the case when the quark multiply 
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scatters from the Color Glass Condensate, then emits a photon and again, mul- 
tiply scatters from the Color Glass Condensate. Those processes are illustrated 
in the figure [l]. Note that Eq. ( [i"3| ) is an exact formula for the bremsstrahlung of 
a quark in a colored glass condensate. Indeed, it resums the interactions to all 
orders with the classical background field, as required by the large gluon density 
in the nucleus. 




Figure 1: The physical processes contributing to the bremsstrahlung of a quark 
in a colored glass condensate. The black dot denotes the interaction of a quark 
to all orders with the background field, i.e. the scattering matrix T F defined in 
Eq. ©. 



Now we will show that the last term also vanishes. To show this, it is 
sufficient to consider the structure of the poles in the propagators in the last 
term. One has the following structure 

+00 

6(q- - r)5(k- + r - p~) ( dl+ J— ,}, . 2 , 2 . (14) 

-oo 1 2F ' +K 2(k-+l-) 

Since both > and q~ > 0, then l~ > as well as k~ + l~ > 0. This 
means that both l + poles are below the real axis. Therefore, one can close the 
integration contour above the real axis and get a vanishing contribution. The 
only exception to this argument would occur if there were terms proportional 
to l + in the numerator which can compensate the l + in the denominator, a 
situation in which the theorem of residues would give a nonzero contribution 
from the half circle at infinity used to close the contour. To investigate this 
case, we write explicitly the terms which contain l + in the numerator of the last 
term. We have 

u{q)l" tit l~u(p) = l + l + u{q)j~j~ i 7^7 _ u(p) . (15) 



Since 7 _ 7~ = 0, equation (15) identically vanishes. Therefore, the last term 
in the amplitude ( |l3| ) is zero. That this term vanishes has a simple physical 
explanation. Indeed, in this term the quark first scatters off the nucleus, then 
propagates for a while and emits the photon, then propagates again and finally 
reinteracts with the nucleus. Note that the intermediate free propagations are 
off-shell and take some nonzero amount of time. But since the nucleus moves at 
the speed of light in the +z direction in our model, it is already far away behind 
the quark by the time the photon has been emitted, and a second scattering of 
the quark off the nucleus cannot happen^} 

d One can notice that the argument used to prove that Eq. (O) vanishes would not be 
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The amplitude can then be written as 
{q(qh(k) ont \q(p)in) = 

T F (q,p-k)Gl( P ^k)^ + ^G F (q + k)T F (q + k,p) 



-eu{q) 

-ieu(q) 
x2n5(q 



u(p) 



(p — k) 2 — m 2 



{q + kf 



p ) J d 2 x±e 



:(q±+k ± -p ± )-x 



u{p) 

*-(u(x ± )-l) .(16) 



Depending on whether one wants to calculate a diffractive or inclusive quantity, 
one would need to color average and then square the amplitude or square the 
amplitude and then color average. In either case, the color averaging does not 
affect the spinor structure and one can evaluate the spinor dependence and color 
averaging independently. 

The spinor dependence of the amplitude (plf) squared is given by 



\ v / /spn 



-tr< (tf+m) 



■ {4- 



j* ^-/t + m ) 7 - ^ 7^(^+^ + m) f 
[p — k) 2 — m 2 



{q + kf 



f + m)$ f (</ + f + m)j 



(p — k) 2 — m 2 (q + k) 2 — 



(17) 



where the factor 1/2 comes from averaging of the spin of the incoming quark. 
Summing also over the spin of the final quark and over the polarization of the 
photon, a straightforward calculation gives: 



(tr(LtL)) 



spin 



-4m 



+8(p~ 



P 



(q-k) 2 (p-k) 2 (p-k)(q-k) 



p ■ q 



1 



(p ■ k)(q ■ k) q-k p-k 



(18) 



3 Color averages 

In order to perform the color averaging of this amplitude with the weight 

* M — >{ - / J } (19, 

we will need to evaluate expressions like (U(x±)) p and (W(x±) U(y±)) where 
U(x±) is given by (|l^). These are already evaluated in |42j1 . Following the 
notations of these papers, we obtain 

(%H), = V(x^)(e-^ 1) (20) 

applicable if the nucleus was moving at a speed v < 1. Indeed, in that case we would not have 
the delta functions implying the conservation of the — component of the quark momentum, 
and we could have poles on both sides of the real i + axis. 



G 



where V(x±) is a function that describes the transverse profile of the nucleus. 
It can be thought of as a function whose value is outside the nucleus and 1 
inside the nucleus. The object B\ appearing in this expression is given by 



Si(xj 



Qt 



d 2 : 



G 2 (x A 



Ql 



A 2 

QCD 



(21) 



with Ql = g 4 (t a t a ) f*™ dz li 2 (z )/2 the saturation scaled] (the integral of ll 2 
over z~ is the number density of color sources per unit of transverse area in the 
target nucleus). Similarly, we have 

((rf(x x ) - l)(U(yj_) - l)) p = P(x ± )V(yx) [l + e~ B ^-^ - 2 e - Bl ] (22) 
with the definition 



B 2 (xj 



Q 2 S J d 2 z ± [G (x 

Ql{x±-y±^ 2 



In 



z±) - G Q (y ± - z ± )] 2 
1 



(23) 



In the above equations, Gq(z± 
defined by 



4tt "\\xj_-y ± \A QC ^ 
y±) is the free propagator in two dimensions, 



d 2 
dz\ 



Gq{z\ 



H z ± - yj 



and given explicitly by 
Gq{z± 



cj_ e' 



■(zjl-Vx) 



(2n) 2 



fc 



(24) 



(25) 



Note that the objects evaluated in Eqs. (gfj) and (|22j) are matrices in the funda- 
mental representation of SU (N c ) that are proportional to the unit matrix. In 
the calculation of cross-sections, one must sum over the color of the outgoing 
quark and average over the color of the incoming quark, which amounts to tak- 
ing the color trace of this matrix and dividing by N c . Therefore, Eqs. (|2^) and 
T2h can be seen as scalars giving directly the result of this procedure. 



4 Cross-Section 

At first sight, the square of the delta function 5(p~ — k~ — q~) that appears 
when we square the amplitude might seem a little worrisome. However, this is 

4 The saturation momentum would acquire a dependence on the rapidity of the quark via 
quantum evolution effects not included explicitly here. Indeed, the quark is sensitive to all 
of the nucleus constituents that have a rapidity between the quark rapidity and the nucleus 
rapidity. 
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just a manifestation of Fermi's Golden rule, adapted to the symmetries of the 
present problem. Indeed, the target nucleus being invariant under (light-cone) 
time £ + , the — component of the projectile momentum is conserved, hence the 
5{p~ — k~ — q~) at the amplitude level. Then, the 2ttS(0~) we have in the 
amplitude squared is just an artifact of shooting a plane wave at this target 
instead of a properly normalized wave packet. It should be interpreted as a 
large but finite time Ax + (see [Q, pages 96 — 97 and |5Cfl, pages 100 — 102 for 
a nice discussion of this). It is in fact simpler to follow pl| (pages 99-107) and 
introduce a wave packet instead of the incoming plane wave 

f d 3 l P ibl ± 

|& n ) EE / ° „ *-- </>(!) g(0in) , (26) 



where 4>(l) is some wave packet that is peaked around the central value p. The 
phase cxp(zb • ij_) is here to account for a finite impact parameter between 
the wave packet of the quark and the trajectory of the center of the nucleus. 
Normalization is chosen in such a way that: 

(bn\<t> ta ) = l ) i.e. I ^L|0(O| 2 = 1. (27) 

Then, the differential (per unit of volume in the invariant phase space of the 
final state particles) interaction probability between this wave packet and the 
nucleus is simply 

dP ' 6 ' 5 (2^(2^J <^)-'l^>| 2 - < 28 » 

Recalling now the relation between the differential cross-section and the proba- 
bility defined in the previous equation 

da = [ d 2 b dP(b) , (29) 



we have for the cross-section: 

d s k d* q r dn d?v ib . (l± _, ±) m 

(2^2fc (27r)32g J {2nf {2nf ^2El ^ 

X(«(<l)7(fc)aut|«(0ta>(«(i')in|«(«l)7(fc)out> • (30) 

Performing the integration over the impact parameter b brings a (2ir) 2 S(l± — 
l'x)- If w e factor out the delta function contained in the amplitude in the 
following way 

(g(gfr(*0aut|g(Oin) = 2*6(r - q~ - k~)M(l\qk) , (31) 

we can also perform the integral over d 3 l' by using one of those delta functions 
and get a factor ^2E p /2p~ , where we have taken advantage of the fact that 
the wave packets are peaked around p to replace by p the dummy variables I 



and I' whenever they are in slowly varying quantities (i.e. everywhere except in 
the wave-packets themselves). At this stage, the differential cross-section reads 

{2ir) i 2k Q (2n) 6 2q 2p 
This is the formula we are going to use in the rest of this paper. 



5 Diffraction 



In order to calculate the diffractive cross section, we average the amplitude ( pL6| ) 
over color charge and then square it. This involves evaluating expressions like 
(U(x±) — 1) which is given by Eq. (20). Performing also the integrals with 
respect to the transverse positions, we obtain 



\M(p\qk)\l 



diff 



(tr(iti)) 



spin 



xV(p±_ -q±- k±)V(~PA 



k 



(33) 



where the function P(l±) is the Fourier transform of the profile function V(x±). 
Note that for a large nucleus, this Fourier transform is very sharply peaked 
around l± = 0, with a typical width of 1/R where R is the radius of the nucleus. 
Since this momentum scale is much smaller than any other typical momentum 
scale in the problem, we can approximate 



V( Pl 



k±)V(-p A 



fei 



V(0)(2tt) 2 6(p a 
ttR 2 (2tt) 2 S(pj_ 



fei 



(34) 



Therefore, we can rewrite the square of the diffractive amplitude as 



\M(p\qk)\ 2 dis = e 2 7rR 2 (tt&L))^ 



1 - e 



(27r) 2 <5(p ± -g ± -fc ± ) , (35) 



so that the diffractive cross-section for bremsstrahlung becomes: 



afJ diff — 



d 3 k 



d 3 q e 2 irR 2 



(27r) 3 2fc (27r) 3 2g 



x(27r) 3 ,5(pj 



k±)S(p~ 



spin 

q--k~) 



(36) 



To proceed further, we assume that the incoming quark is massless and that 
its transverse momentum p± is zero. This is reasonable if we have in mind an 
application of this calculation to photon production in pA collisions. Indeed, 
there is no significant amount of heavy quarks in the wave function of the proton 
and partons are collinear to the proton. Under these assumptions, it is a matter 
of straightforward algebra to show that 



(tr(LtL)) diff = 

\ v } I spin 



(37) 
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leading to 



%S--0. (38) 



This result is in fact rather intuitive: asking for a diffractive process implies 
that no net transverse momentum is exchanged between the classical color field 
and the quark, which inhibits the radiation of a photon by bremsstrahlung. 



6 Inclusive 

In order to calculate the inclusive cross section, we square the amplitude and 
then average over color charge. Using (|2^), we obtain for the inclusive amplitude 
squared 

\M(p\qk)t d = e 2 {tr(L^L)) sp . in Jd 2 x^d 2 y^e^ +k --^y^-y^ 

xV(xA_)V{y±) \l + e- B ^~ y ^ - 2e~ m ] . (39) 

By inspecting the formula (|23|), we see that exp(— B^) becomes very small if 
\x± — y±\ becomes larger than 1/Q S . Therefore, this term contributes only 
if the separation between x± and y^ is much smaller than the radius of the 
nucleus (assuming a large nucleus such that R ^> Qj 1 )- Therefore, for this 
term we can approximate 

P{x±)P(y ± ) « P 2 (x±) « P{x±) , (40) 

where the last equality is valid for a nucleus with rather sharp edges (i.e. if 
V(xj_) goes from to 1 rather fast). Using Eq. ( p4] ) and performing the Fourier 
transforms, we can rewrite the inclusive amplitude squared as follows 

|M(p| 9 fe)lLl = eW(tr(it L )) sp . n 

x [D(p ± -q ± - fe x ) + (1 - 2e- B ')(2n) 2 6( Pl _ -q ± - k±)] , 

(41) 

where, following p^| , we denote 

D(l±) ee J d 2 x ± e ll ^ x ^{e~ B ^ x ^ - 1) = J d 2 x ± e Ux3!x (U(Q)U*(x±) - l) p . 

. (42) 

Inserting this in the formula (|3^) , we obtain the following formula for the inclu- 
sive cross-section 5 

5 The invariant phase space of a massless particle can be rewritten in terms of the light-cone 
variables as follows: 

d 3 k dk+dk~d 2 k± „+.,., j_ 



(27r) 3 2fc (2tt) 4 



27r6»(fc ± )<5(2fc-fc+ - k\) . (43) 
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x [D(p ± -q ± - fe x ) + (1 - 2e- B ^(2n) 2 S{ Pl _ - q x - k±)] . 

(44) 

Another simplification can be achieved by neglecting the term in exp(— B\) since 
Bi ~ Q 2 /K 2 D 3> 1 appears in the exponential with a negative sign. If one 
introduces | p5[ 

C(I X ) = [ d 2 x x e il ^e- B2 ^ = [ d 2 x ± e il ± *± (U(0)U\x ± )) , (45) 



the inclusive cross section can be rewritten as 



dCTi " cl _ (27r)3 2 fc (271)32^ 2p- 1 ( ^ 
x27r5(p _ — g~ - fc~) C(px - <Zx — fei.) • 



(46) 



Assuming again that the incoming quark transverse momentum p± is zero and 
neglecting the quark mass, one can perform the integrals over q + ,k + , q~ using 
the delta functions. There is however a complication due to collinear singular- 
ities, i.e. singularities that show up when the emitted photon is parallel to the 
outgoing quark. It is convenient to trade the transverse momentum of the final 
quark for the total transverse momentum of the final state, i.e. l± = q± + k±. 
In terms of this new variable, we have 

= 2 r toZ+SL^l [ A± JlcM ,47, 



nR 2 d 2 kj_ (27r) 5 fc;[ J z J [lx-kx/z] 

where z = k~/p~ and [1 + (1 — z) 2 ]/z is the standard leading order photon 
splitting function. Eq. fl47| ) is our main result. Note that C(l±) behaves like 
at large l± which ensures that the integral converges at large momentum 
transfer. In this formula, C(l±) is the only object that depends on the struc- 
ture of the color sources describing the target nucleus. In particular, all the 
quantum evolution effects would go into this object via the averaging procedure 
in Eq. (pOf). One can also note that this result exhibits the standard collinear 
denominator [l± — k±/z] 2 that vanishes if the photon is emitted collinearly to 
the quark. This aspect of the result is of course not affected by the description 
of the target nucleus as a color glass condensate. 

In the soft photon limit, one can see the decoupling of the photon emission 
subprocess from the quark scattering part. The latter agrees with the quark- 
nucleus scattering cross-section calculated in Q . 

It is instructive to perform the "perturbative limit" of this result. This 
regime is reached when the transverse momentum l± transferred between the 
nucleus and the quark is large compared to the saturation momentum Q s . In 
this limit, we have Ell 

20 2 

C(l±) « • (48) 
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Figure 2: Behavior of the correlator C(k±) as a function of k±. In this plot, 
the value of Q s is such that Q S /A QCD = 10. Circles: computed value of C(k±). 
Solid line: the "perturbative" value C(k±) « 2Q s /k 4 L , valid if k± 3> Q s . 



Using this result, we have: 



da; 



incl 



1 



d 2 k 1 _ 



2N h e 2 q a em a 2 C F f 1 + (1 - z) 2 f d 2 l 



pert. IT 2 k\J ^ Z J Il[l±~k ± /Z} 2 ^ 





where e q is the quark electric charge in units of the electron charge, and where 
Nh = ttR 2 J dz~ fj 2 (z~) is the total number of hard color sources in the target 
nucleus. Therefore, this expression has all the features of the bremsstrahlung of 
a photon by a quark scattering off a parton inside the nucleus with the exchange 
of a gluon in the t-channel (this term is the dominant one at large center of mass 
energy). 

In Eq. (p7|), the only factor that depends crucially on the saturation hypoth- 
esis for the nucleus is the factor C(l±). Indeed, this term contains all the de- 
pendence on the saturation scale, as well as the modifications of the transverse 
momentum spectrum at scales below Q s . The transverse momentum depen- 
dence of this object is illustrated in figure 0. In order to observe effects due to 
this factor, it would be useful to measure both the radiated photon and the jet 
induced by the outgoing quark. The photon-jet correlations, and in particular 
the distribution of their total transverse momentum l± — q±+kj_, would indeed 
enable one to extract in a rather direct way the function C(l±) itself. On the 
contrary, if one measures only the photon spectrum, one can access only a given 
moment of this function. 
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7 Conclusions and perspectives 



In this paper, we have studied the photon production in pA collisions by bremss- 
trahlung of a quark in the saturated color field of the nucleus. We have shown 
that the cross-section for this process depends sensitively on saturation effects. 
Furthermore, we have shown that photon bremsstrahlung diagrams which are, 
in the standard perturbation theory, higher order compared to the direct photon 
diagrams, become leading order. This is due to the presence of a strong color 
field 0(1/(7) in the target nucleus which compensates for the factor of g from 
quark-quark-gluon vertex. 

Experimentally, it may be an easier task to look for a photon and a jet pair 
rather than a single jet. This was studied in the context of heavy ion collisions 
in p2j. Therefore one can directly measure the correlator C(l±) that describes 
the interactions of a high energy probe with the target nucleus. It should be 
noted that saturation effects become larger in the forward region and therefore 
easier to measure experimentally. It is straightforward to extend this calculation 
to the case of production of dileptons and di-photons. This is currently under 
investigation. 

Another interesting extension of this work would be to study the radiation 
of a gluon by a quark as it moves through the classical field of the nucleus. 
In that case, the emitted gluon can also interact with the background classical 
field, which makes the calculation of this process more involved. 
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